Solid state effects on exciton states and optical properties of PPV 
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We perform ab initio calculations of optical properties for a typical semiconductor conjugated 
polymer, poly-para-phenylenevinylene, in both isolated chain and crystalline packing. In order to 
obtain results for excitonic energies and real-space wavefunctions we explicitly include electron-hole 
interaction within the density-matrix formalism. We find that the details of crystalline arrangement 
crucially affect the optical properties, leading to a richer exciton structure and opening non-radiative 
decay channels. This has implications for the optical activity and optoelectronic applications of 
polymer films. 
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Ordered films of organic conjugated polymers are of 
strategic relevance for novel optoelectronic devices [Q. 
In addition, they offer an ideal scenario for the study of 
electronic and excitonic confinement, being composed of 
quasi-one-dimensional (ID) systems, arranged however in 
a three-dimensional (3D) crystalline environment. Long 
linear chains can indeed be thought of as ID systems 
with the highest degree of ID confinement for a poly- 
atomic system (~ 2-5 A). Since the main optoelectronic 
characteristics derive from the mobile 7r-electrons, delo- 
calized along the chain backbone, and non bonding to 
neighbouring chains, the vast majority of studies adopt 
the single-chain model: in fact, complete ab initio the- 
oretical studies of the optical properties have been per- 
formed for isolated chains, highlighting the strong con- 
finement expected for such systems H . Recently, very 
simplified models of polymers in a "crystalline medium" 
have been approached through ab initio theory Q ; how- 
ever the effect of crystal structure or side chains is com- 
pletely neglected. This is also usual practice in the in- 
tepretation of experimental data: e.g. data on poly-para- 
phenylenevinylene (PPV) and its alkylated and metoxi- 
lated derivatives (MH-PPV, MEH-PPV, etc. . . ) are usu- 
ally bundled together as representative of PPV, based on 
the supposition that details of the 3D structure are not 
relevant. The picture emerging from this analysis is far 
from clear, however: quoted exciton binding energies dif- 
fer by an order of mag nitude |, |, |, |, |, [0|, there 
is an on-going controversy about the existence of charge- 
transfer excitons or excimers ||, ^ |l^, |l4j , and about 
conditions for efficient light emission. The question then 
arises if, on the contrary, solid state effects cannot be ne- 
gleted and crystallization entails also structure-specific 
3D interchain coupling effects: recent data on oligoth- 
iophene crystals seem to suggest that this might be the 
case H 

Here we address these issues through a full ab initio cal- 
culation of optical spectra and real-space exciton wave- 
functions for PPV, in both the isolated chain and the 



crystalline phase. In the single PPV chain, we find that 
the lowest singlet exciton extends over a few monomers 
along the chain and is optically active, with a binding 
energy of about 0.6-^-0.7 eV. In the crystal we find ev- 
idence of significant interchain coupling, with decrease 
of the first optically active exciton binding energy to 
about 0.2-^0.3 eV |0, ^, ^; as a consequence of both 
interchain interaction and detailed crystalline structure, 
a dark split-off singlet exciton appears at lower energy, 
which introduces a new path for non-radiative emission 
and electron-hole dissociation, and is expected to de- 
crease drastically the photoluminescence efficiency [ p^ ; 
furthermore, we find a much richer exciton structure, 
with charge-transfer excitons close enough to the first 
direct exciton Q to account for early electron-hole dis- 
sociation and apparently- vanishing binding energies [p^ ; 
fitting all the pieces together, our results can explain the 
controversies on the optical behavior of these crystals. 

The exciton binding energy, Eb, is a central quantity 
in the photophysics of polymer semiconductors because 
it is related to the spatial extension of the excitation and 
to the probability of radiative emission/ absorption and 
electric-field-induced dissociation of the exciton. It is de- 
fined as the energy required to separate a bound electron- 
hole pair, the exciton, into a free electron and a free hole 
0, |, H |ll|, [l^- For PPV, the experimental values that 
are currently reported range from 0.2-^0.4 eV @, ||, ||] to 
around 1 eV [§Jq|, but energies of less than 0.1 eV are 
also proposed |10(]. These values correspond to spatial 
extensions varying from a single monomer (Frenkel-like 
excitons) up to many unit cells (Wannier-like excitons). 

From the theoretical point of view, calculating the 
properties of excitons in extended polymeric systems 
is much more difficult than in conventional inorganic- 
semiconductor quantum wires |]l8j, because it requires 
to combine an accurate treatment of the Coulomb in- 
teraction with a realistic microscopic description of the 
single-particle electronic structure beyond effective mass 
approximations. Standard quantum-chemistry schemes 
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|l9[ for treating excitations in finite molecular systems 
are intrinsically not suitable for crystals. To this end, 
we have implemented a first-principles density-matrix ap- 
proach where the correlated electron-hole spectrum is ob- 
tained directly from the interband polarizations {d^cp, 
resulting from the creation of an electron in conduction 
state /i, with energy and a hole in valence state v with 
energy ejj (cj^ and dj, are the Fermionic creation opera- 
tors; the labels ^ = (c, k) and v = (w,k) stand for both 
the band index and the wavevector). The polarization 
eigenstates are obtained by direct diagonalization of the 
two-body Schrodinger equation: 
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Here the first term on the right-hand side accounts for 
the uncorrelated electron and hole quasi-particle ener- 
gies, the second for electron-hole screened direct (W) and 
unscreened exchange {V) Coulomb interaction {Ss = 1 
for singlets, otherwise). The interaction kernels are 
written in terms of the single-particle wavefunctions ip: 

W^MM>.' = / {viW^. (ri )u;(ri , rs) V^* (rs)^^ {v2)dv^dv2 



r:{vi)i^t.{ri)v{r,,r2Wi*{v2W^.{v2) dv^dv^ 



where u(ri,r2) — \r\ — is the bare Coulomb po- 

tential while w(ri,r2) — ^ e~^{y\,v)v{y,v-i)dv is the 
screened one. The excitonic eigenenergies Ex and eigen- 
states A are provided by the solution of the above eigen- 
value problem; hence, this description closely re- 
sembles the earlier one obtained within a Green's func- 
tion scheme, which comprises the solution of the Bethe- 
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therefore the calculation of non- 



Salpeter equation |21 
The first 
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interacting electron and hole quasiparticle energies e 
and wavefunctions tl) that we obtain through a density- 
functional theory hamiltonian in the local density ap- 
proximation [ p3| , p4[ (with norm- conserving pseudopoten- 
tials and a plane- wave basis with 50 Ry energy cutoff), 
and by subsequently adding a rigid shift to the conduc- 
tion band energies to account for the self-energy correc- 
tion (taken from Ref. ||). Moreover, in order to make 
the full calculation computationally feasible for complex 
systems as polymer crystals, we approximate the full ef- 
fective screening function by a constant diagonal dielec- 
tric tensor e, whose components are calculated ah initio 
within linear response |2^. For the case of the isolated 
polymer chain, the assumption of space homogeneity is 
indeed very drastic, so we introduce an effective screen- 
ing tensor, which is defined within the volume r^eff where 
the interaction actually takes place, while e = 1 out- 
side Oeff. The components of e are again calculated ah 
initio within linear response and fleff is determined self- 
consistently with e to encompass the region of significant 
electron-hole correlation function. 
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FIG. 1: z-polarized optical spectra for single chain PPV 
(left) and the PPV crystal in herring-bone structure (right). 
The solid lines are the 2-polarized spectra calculated including 
electron-hole Coulomb interaction; the single particle spectra 
are shown for comparison as shaded backgrounds. An arti- 
ficial broadening of 0.04 eV (0.08 eV) was included for plot- 
ting the Coulomb-correlated (single particle) spectra. The 
exciton binding energy is the energy difference between the 
exciton absorption peak in the correlated spectrum and the 
onset of the continuum in the single-particle spectrum. The 
arrows mark the lowest singlet exciton in the isolated chain 
spectrum, and relevant singlet excited states in the crystal 
spectrum: first optically active direct exciton (DA), optically 
inactive (dark) direct exciton (DD); charge-transfer excitons 
(CTl and CT2). 



In order to better understand the role of interchain 
coupling for the crystal case, we first consider the iso- 
lated chain as reference system. We show in Fig. |^(a) 
the results for the absorption spectrum of the isolated 
chain j^, which is highly polarized along the chain 
(z-direction); for x and y polarizations absorption close 
to band gap is negligible. We find an optically active 
singlet exciton with energy approximately 0.6-^0.7 
eV below the single-particle band-gap energy. This low- 
est singlet state extends over a few unit cells along the 
chain ^, the large binding energy arising essentially 
from the strong lateral confinement. To gain insight in 
the nature of this state, in Fig. ^ we plot its electron and 
hole contributions: p^{vh) = Yl,yi\Yl,v '^i^v'^^t{'^h)f and 
P''(re) = E,. I ^/^i'V';:(i"e)P- Both have their maxima 
on the vinylene chains but in different bonding config- 
urations, and thus vinylene atoms are expected to un- 
dergo the largest photoinduced lattice relaxation; this 
suggests the assignment of the phonon replica observed 
in PPV photoluminescence to the vinylene C-C stretch 
mode (around 0.2 eV |25|, [2^ , p7| , p8[ ) and explains the dif- 
ference in lineshape |29[] with poly-para-phenylcnc (PPP), 
where the vinylene chains are not present. 

We point out that the exciton signature is accompa- 
nied by suppression of the typical ID van Hove singu- 
larity (^ I/a/E) [§■ This is a phenomenon resulting 
from dimensionality, also found in inorganic semiconduc- 
tor quantum wires It is relevant for the ex- 
perimental determination of the exciton binding energy: 
the energy separation between the exciton and the on- 
set of the ID absorption continuum cannot be extracted 
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FIG. 2: Electron and hole contributions, p'^ (a) and p** (b), to 
the lowest singlet exciton of the isolated PPV chain, plotted in 
a plane parallel to the phenyl rings. The colour code is shown 
on the side bar, with p'^ {p^) normalized to its maximum 
value. 



directly from a one-photon optical spectrum, since the 
oscillator strength associated to the latter is vanishingly 
small [18 ; a separate measurement is required to extract 
the quasiparticle gap H, 0, ||, |3l) . 

We now turn to the PPV crystal, that we calculate 
in a realistic herring-bone structure |32[ |3^, by keep- 
ing an orthorombic cell with a — 8.07 K, h — 5.08 A, 
c — 6.54 A and the setting angle $ = 52° (see inset 
in Fig. H). The influence of crystalline structure can be 
felt already on the single-particle bands: due to the spe- 
cific crystal symmetry, with two chains per unit cell, all 
bands are doubled |3^; the splitting between doublets 
gives rise to the splitting between excitonic states (the 
so-called Davydov splitting [Q), and is strongly depen- 
dent on interchain interaction. This feature alone points 
to a clear difference between PPV and other substituted 
PPV derivatives, such as MEH-PPV, that cannot pack 
in the herring-bone structure, have only one chain per 
cell, and thus would not show Davydov splittings. Note 
that this very important feature cannot be seen either 
from structureless crystal models |^ or simple molecular 
aggregate models without translational symmetry [^ . 
Other general features are that dispersion along the F-Z 
direction (along the chain) is of the order of ~ 2.2 eV, 
as for the isolated chain, and dispersion perpendicular to 
the chain is vanishingly small. 

We show in Fig. |^(b) the correlated spectrum for the 
crystal. It is again characterized by a strong singlet ex- 
citon peak (marked DA in the figure) and quenching of 
the oscillator strength at the onset of the single-particle 
continuum. While the single-particle band gap is reduced 
in the crystal |Q , the optical gap in the correlated spec- 
trum is close to the single-chain result. Hence, the bind- 
ing energy of this exciton is reduced to about 0.2 eV, 
consistent with the increased screening in the crystal . 
The first absorption peak is still strongly polarized along 
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FIG. 3: Probability distribution \^'='"'{r^,rh)\^ of the lowest 
singlet excitons in the PPV crystal. The state (a) is the op- 
tically active direct exciton DA appearing in the spectra of 
Fig. hi, while the state (b) corresponds to the dark state CTl, 
at slightly higher energy. Here VJ(<==='= is plotted as a function of 
the electron-hole distance in the xz plane (red line in the inset) 
averaged over the center-of-mass coordinate R = (re +rh)/2. 
The colour code is shown below the figure, with the probabil- 
ity normalized to its maximum value in each panel. The inset 
shows a representation of the chains projection on the (a,b) 
plane of the orthorombic cell. 



z, confirming the strong optical anisotropy of the mate- 
rial. 

Now we turn our attention to the other remarkable 
features of the optical behavior: (a) the single-particle 
band-splitting results in an exciton Davydov splitting of 
~ 120 meV, with the corresponding introduction of an 
optically inactive (dark) state below the first active exci- 
ton (DD arrow in the figure); (b) other excitonic states 
are also introduced above the DA peak, with no oscillator 
strength along 2, but just barely active on the perpendic- 
ular polarization (CTl and CT2 in the figure). 

The corresponding exciton wavefunctions for DA and 
CTl are shown in Fig. |^, where we plot the probabil- 
ity density distributions \'^'^^'^{re,rh)\'^ ■ Note that the 
exciton wavefunction ^*^^'^(re, r/i), expressing the prob- 
ability amplitude of finding the electron in r^. and the 
hole in Vh, is computed within our formalism directly in 
terms of the interband polarization A as \l/'^^'^(re, r^) — 
I]^i.^Mi''^Jl(re)V'J^*(r/i)- The DA state is remarkably 
confined on a single chain, only slightly more extended in 
z-direction than in the isolated chain. (The DD state has 
the same distribution as the DA, by symmetry: they are 
both direct as electrons and holes are on the same chain.) 
The CTl state, on the other hand, shows very small on- 
chain contribution, with high probability of finding elec- 
tron and hole on different chains, and can be character- 
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ized as an excimer; same comments apply to the CT2 
state. 

In summary, we have found that the lowest correlated 
states of PP V crystal do involve more than one chain and 
are very dependent on crystal structure [^; however, 
they retain a quasi-one-dimensional character, and the 
consequent suppression of ID singularity precludes the 
use of one-photon experiments to extract exciton bind- 
ing energies. Furthermore, we have found that, due to 
interchain coupling and symmetry effects, the lowest elec- 
tronic excitation is a dark exciton, that is optically inac- 
tive but provides a non-radiative decay path. This can 
strongly quench the photoluminescence with respect to 
non-interacting chains. In view of the impact of crys- 
tallization on the optical properties, we studied also the 
influence of the packing density: without changing the 
symmetry, we increased the interchain distance in the 
TT-stack direction b (see inset in Fig. we find that an 
increase of ~ 20 % leads to an increase in the binding en- 
ergy (to ~ 0.4 eV) and a drastic decrease in the Davydov 
splitting by an order of magnitude, thus restoring ra- 
diative efficiency. These results suggest that the spread 
in photoluminescence (PL) and photoconductivity (PC) 
data observed by different groups may arise from differ- 
ences in the aggregation state of the polymer. Moreover, 
since real films have regions with different chain pack- 
ings, the dominant contributions to PL and PC can orig- 
inate from different regions of the same sample (the least 
and the most closely packed, respectively). Another solid 
state effect is the appearance of bound charge-transfer 
excitations close above the first exciton peak: this not 
only can enhance photoconductivity, but also can be the 
origin of the infra-red active modes seen on two-photon 
^ B6|. 
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experiments 

Our results show that solid state effects are indeed very 
important for the final optical structure of a polymer 
film. This conclusion could only be reached because of 
the possibility of calculating first-principles optical spec- 
tra of realistic three-dimensional polymer structures. 
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